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1 . Introduction 


The work being performed under this grant Is directed towards the 
Improvement of our knowledge of tiu gravity field of the earth and In 
the analysis of quantities that depend on the gravity field. 


2. Work Accomplished 

Under this grant research has been conducted In several different 
topic areas. This research has been carried out by Graduate Research 
Associates and the Principal Investigator. We now describe the activities 
that have taken place during this reporting period. 


2.1 Terrestrial Anomaly Field Improvanent 

We have continued during this period to update our 1^^ x 1<^ terres- 
trial me. n anomaly tape. The version developed during this reporting 
period Is called the October 79 version. This was an update of a June 
78 data tape. The new tape contains 41973 anomalies based on data from 
45 data sources. Of these anomalies 721 values had standard deviations 
greater than 30 mgals and 16209 values were based on data sources exclusive 
of the Defense Mapping Agency Aerospace Center. In addition we Identified 
6912 anomalies that may have been predicted through geophysical correla- 
tion techniques. Additional details may be found In Rapp (1980a). 

Frcmi this data we have also prepared two global sets of 1654 5 ° 
equal area anomalies. The first set was based on all 1° x 1^ data exclud- 
ing anomalies whose standard deviations were greater than 30 mgals. The 
second set was similar to the above but deleted the geophysically predicted 
anomalies . 

The 1° X 10 and 50 anomaly sets were sent to Frank Lerch on January 9, 
1980. We also sent the combined terrestrlal/altlmeter data file con- 
taining a merger of the October 79 tape with the anomalies devised from 
the Geos-3 altimeter data. 

We have also prepared a document that explains our updating procedures 
and the computer programs Involved. The document Is entitled: Programs 

and Procedures for Updating the OSU 1° x 1® Gravity Anomaly Tapes. A copy 
of this report was sent to Frank Lerch on March 5, 1980. 

Currently we are evaluating several new data sources. We tried at 
the lUGG/IAG meeting in Canbera to establish some data exchange with 
some Chinese colleagues but to this point nothing specific has happened. 


2.2 Improved Cfeold Determinations 

Research was also carried out In seeking better analytical techniques 
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for geoid undulation cooputatlons . The procedure to be iiaproved upon 
is the case where potential coefficients and surface gravity anomaly 
Information are combined to determine detailed geoid information. 

The results of this research have been quite successful and are 
documented In three places (Rapp, 1979, 1980b) and Jekell (1980). The 
abstract of the Jekell report Is as follows: 

The truncation theory as It pertains to the calculation of geoid 
undulations based on Stokes' Integral, but from limited gravity data. 

Is reexamined. Specifically, the improved procedures of Molodenskil 
et al. (1962) are shown through numerical Investigations to yield sub- 
stantially smaller errors than the conventional method that Is often 
applied In practice. In this improved method, as well as In a simpler 
alternative (Melssl, 1971b) to the conventional approach, Che Stokes' 
kernel Is suitably modified In order to accelerate the rate of conver- 
gence of the error series. These modified methods, however, effect a 
reduction In the error only If a set of low-degree potential harmonic 
coefficients Is utilized In the computation. Consider, for example, the 
situation In which gravity anomalies are given In a cap of radius 10° 
and the GEM 9 (20,20) potential field Is used. Then C 3 mlcally, the 
error In Che ccxBputed undulation (aside fr<»n Che spherical approximation 
and errors In the gravity anomaly data) according Co the conventional 
truncation theory Is 1.09 m; with Melssl' s modification it reduces to 
0.41 m, while Molodenskil 's Improved method gives 0.46 m. A further 
alteration of Molodenskil ' s method Is developed and yields an RMS error 
of 0.33 m. These values reflect the effect of the truncation, as well 
as the errors In the GEM 9 harmonic coefficients. The considerable 
improvement, suggested by these results, of the modified methods over 
the conventional procedure Is verified with actual gravity anomaly data 
In two oceanic regions, where the GEOS-3 altimeter geoid serves as the 
basis for comparison. The optimal method of truncation. Investigated 
by Colombo (1977), Is extremely Ill-conditioned. It is shown that with 
no corresponding regularization, this procedure Is Inapplicable. 


2.3 GRAVSAT Analysis 

Although not a formal part of our grant effort, some research was 
carried out In an error analysis approach to considering the accuracy of 
the determination of mean anomalies and mean undulations from a GRAVSAT 
mission. The analysis carried out was based on considering error propa- 
gation errors, and truncation errors. We now have a theory and an 
operating program to handle many different mission parameters. Some 
results are described In Rapp (1980a) . A detailed report by Jekell and 
Rapp Is now being prepared. 
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3 . Future Work 


We will coQCiaue our 1° x 1® data collection work. More Intensive 
work is needed to identify bad 1® x 1® anomalies. For our next update 
we. ma' >.<re additional Ceos-3 and Seasat data that may be used to verify 
grouni ' ased data. 

We will continue to work with our improved geoid computation proce- 
dures with specific attention being given to using high degree reference 
fields such as GEMIOB. The problem here is to obtain a realistic accuracy 
assessment of the potential coefficients. 


4 . Personnel 


During this reporting period the following persons have narticipated 
in the work of this grant: 

Richard H. Rapp, Project Supervisor 
Thomas Croxell, Graduate Research Associate 
Christopher Jekeli, Graduate Research Associate 
Kostas Katsambalos, Graduate Research Associate 
Agit Singh, Graduate Research Associate 
Kamil Eren, Graduate Research Associate 
Pamela Pozderac, Secretary 
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We give here reports and papers that have been produced under this 

grant since the last status report. 

Christodoulidis , D., Accuracy of Relative Oceanic Geoid Computations, 

J. Geophy. Res., 84, Sept. 10, 1979. 

Christodoulidis, D., Influence of the Atmosphere Masses on the Gravita- 
tional Masses of the Earth, Bulletin Geodeslque, Vol. 53, No. 1, 
1979. 

Jekeli, C. , Accuracy Estimates of Geoid Undulation Differences, Computed 
from Gravity Disturbances, iVnomalies and Potential Coefficients, 
paper presented at the 1979 Spring Meeting of the American Geophysi- 
cal Union, Washington, D.C., May 1979. 

Jekeli, C., Global .Accuracy Estimates of °oint and Mean Undulation Dif- 
ferences Obtained from Gravity Disturbances, Gravity Anomalies and 
Potential Coefficients, Department of ceodetic Science Report No. 
288, May 1979. 

- - - - Programs and Procedures for Updating the OSU 1® x 1® Gravity 
Anomaly Tapes, Internal report. Department of Geodetic Science, The 
Ohio State University, February, 1980. 
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Katsa^Ios, R. , TI'C Effect of the Saoochiog Operator on Potential 
Coefficient Detendnations . Departa«it of Ceodecic Science, 

Report No. 287, March 1979. 

Rapp, R. H. , The Oceanic Ceoid, paper presented at the X7II General 
AsscBbly, IUGG,IAG, Canbera, Australia, 1979. 

Rapp, R. H. , The 1<> x 1° Hem Anoaaly Field of tiw Earth and Prospects 
for Its laproveaent, paper presented at the 4th Intemational 
S 3 ;iBposiua, Geodesy and Physics of the Earth, Rarl-Marx-St»lt, 

GDR, May 1980. 

Rapp, R. H. , A Comparison of Altiaieter and Gravimetric Geoids in the 
Tonga Trench and Indian Ocean Areas, Bulletin Geodesioue, Vol. 

54, 1980b. 


